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Abstract 

Travelling wave in a helical wave guide is considered for electron 
acceleration. A first determination of the travelling wave modes using 
a partial wave expansion (PWE) and a point matching method is 
presented. It gives a rapid solution for moderate deformation of the 
guide relative to a straight cylinder. Strong deformations will give 
higher accelerating gradient but the PWE diverges. Two methods 
overcoming this difficulties are suggested. 



1 Introduction 

An electromagnetic wave in a straight wave guide cannot accelerate particles 
because its phase velocity f ph is larger than c. In conventional linear accelera- 
tors, f ph is reduced by a chain of resonant cavities along the guide. A scheme 
not using a resonance phenomenon is the travelling wave tube (TWT). One 
type of TWT is a cylindrical tube with a conducting helical wire inside [H [2] . 
Another type is a tube with a thick dielectric wall P]. 

In this paper we study another kind of helical device : the helically bent 
tube. The main "photon path" is the helix passing trough the centres of the 
local cross sections of the tube. For a strong enough bending, the difference 
in length between this path and a straight path is such that the effective Uph 
is smaller than c, although the local Vph is larger than c. 

One may expect several advantages in using helically bent tube in place 
of accelerating cavities : easy construction, easy cooling, easy maintainance 
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of the wall surface using an helical arm, large frequency spacing of the modes, 
large group velocity which allow a fast filling of the guide by the wave. Al- 
ternatively, the wave can be generated by a driving beam. Like with the 
dielectric wall in [1] , several modes can cooperate to transfer the energy from 
one beam to the other. 

The helical wave guide is defined by its axis Oz, its period A = 27r/g 
and the base curve, which is the section of the wall with the 2; = plane, 
of parametric equation rfe(s) = {xb{s),yh{s)). The point (x,?/) = (0,0) is 
supposed to be inside the base. The wall surface is generated by helices Hs 
of parametric equation 



where the 2x2 matrix TV^qz] describes the rotation of angle qz about O2;. One 
method to make a helical waveguide is to draw helical grooves on the internal 
face of an initially cylindral tube. The case of rectangular grooves has been 
studied in [5l |6] . The corresponding wall base curves have indentations. Here 
we will restrict to base curves with smooth curvatures. 

This work is the continuation of a first study [7J, where the electromag- 
netic wave modes were calculated using a partial wave expansion (PWE) 
associated to a point matching (PM) method. Since that time, we realized 
that PWE is limited to weak enough helical bending and sufficiently "round" 
base curves ; if it is not the case, the PWE converges only inside a cylinder 
which does not contain the whole guide. 

In this paper we will recall the PWE-PM method, some of our old results, 
present some new ones and discuss the limitations of this method. Two 
alternative methods will be evoked. 

2 Propagating modes in the scalar field case 

We first present the simpler case of a scalar field $(/!:, x, y, 2;), which obeys 
the Klein-Gordon equation (in units c = 1) 



with the boundary condition (^{t,x,y,z) = on the wall. These two equa- 
tions are invariant under a time translation and under a helical displacement 



{x,y,z) = (n[qz]rbis), z) 



(1) 



{-dl + dl + dl + dl) $(t, X, ^) = , 



(2) 



{translation by A2;} x {rotation by qAz} . 



(3) 
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In the transverse plane we introduce a "rotating frame", the coordinates 
(X, Y) of which are related to the fixed frame ones hj x + iy = e*^^ (X + iY). 
On can therefore look for a solution of the form 

<l>(t, X, y, z) = e^^^-*^* ^(X, Y) (4) 

which has definite reduced frequencies u and definite eigenvalue P of the 
operator 

^hoiic =Pz + qJz = -id/d^ - iqd/d^ . (5) 

2.1 The partial wave expansion (scalar field case) 

Let us temporarily assume that a mode can be expanded in free cylindrical 
partial waves which are eigenstates of and Pz with eigenvalues / and p 
respectively. A single partial wave writes 

<l>^,p,Kt, X, y, z) = e-'^' e^P^ e*''^ V'/.fclr) , (6) 

where (r, (^9, 2;) are the cylindrical coordinates and 

^lA"r) = ■h\{^r)|k\'\=^\{Kr)|^^'\. (7) 

Ji and li are the regular Bessel functions of first and second kind, k = 
(cu^ — + iO)^/^ is the radial momentum and n = —ik = [p^ — up' — iO)^/^. 
The two expressions of ([7]) are equivalent. If one whishes to work with real 
quantities, one takes the first one for \p\ < \u\, the second one for \p\ > \u\. 
The denominator fc''' is introduced to make ^i,k{i^) finite at /c — )■ 0. The 
partial wave expansion reads 

+00 

<l>(t, X, y, z) = e-'^' Yl ^'"^'^'^''^ ^'''^(O(^) ' (8) 

l=—oo 

with 

p{l) 01 pi = P-lq, k{l) or ki = {u^ -pf + ^0) , (9) 

so that each partial wave is an eigenstate of "Pheiic with eigenvalue P. 

It is sufficient to write the boundary condition at t = z = 0, thanks to 
the invariances under time translation and helical displacement. For every 
point rb{s), of polar coordinates {rb{s) , ipb{s)) , of the base curve one has the 
constraint 

ai exp[ilipb{s)] V'i,fc[rfe(s)] = . (10) 
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Figure 1: Matching points on the base curve. 



2.2 The point matching method in the scalar field case 

One keeps only a finite set L of partial waves and selects an equal number of 
matching points ri, ...r^v of the base curve (see Fig J2.2l) . For a base without 
discrete rotational symmetry, £ = [/inf , /sup] with Zsup = ^inf + — 1. The 
choice of r„'s and L should be guided by the following conditions : 

• the r„'s should give a precise enough definition of the base curve, 

• £ should contain /'s of both signs, however with more Ts of the same 
sign as q. 

• / (v^n+i — V^n) for / G £ should not be too large compared to unity, 

• ki (r„+i — r„) for / G £ should not be too large compared to unity. 
The ai are approximatively determined by the system of linear equations 



(r„, i^n) being the polar coodinates of r^. A nonzero solution exists if the 
determinant est zero: 





(11) 



det{C„i} = , where Cni 



(12) 
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This equation, which relates u and P, is the dispersion relation in the PWE- 
PM approximation. One may impose u and look for the roots of fll2l) as 
a function of P, or vice-versa. Alternatively, in view of accelerating highly 
relativistic electrons one may impose u/P = fph = f e — 1 and look for the 
roots of f|T2|) as a function of u. Using a graphical software, the roots are 
easily seen as dips in the curve of In |det{Cn;}| versus the free variable. If 
(det{C„/}) is e*" times a real number, the roots are precisely located as jumps 
of the curve arg (ie~*" det{C„/}) from ±7r/2 to =F7r/2. 

Provided that the partial wave expansion converges, the quality of the 
approximation is improved by increasing A^. The convergence is obtained 
first for the lowest modes (which have smallest uj at fixed P or fph, or largest 
P at fixed w), then for higher and higher modes. A first criterium is the 
stability of the root when A^ is increased by one or two units, when £ is 
translated by one or two units or when one moves the matching points. The 
ultimate criterium is the stability of the coefficients a^, at least for I not close 

to /inf or /sup- 

Case of a symmetric base. If the base is symmetrical about the x-axis, 
one can choose the {r^} to be symmetrical by pairs, except for those which 
are on the x-axis. Then i"dei{Cni} is real, v being the number of pairs. The 
ai can be taken all real, which gives 

$(0,x,-y,0) = <I>*(0,x,i/,0). (13) 

If the base is symmetrical about Oz, the modes are either even or odd 
under the parity symmetry {x,y) — > (— x,— y), therefore £ has to contain 
either even or odd /'s. More generally, if the base has a rotational symmetry 
of order S*, the Ts should be equal modulo 5*. To avoid redundant equations 
in ( !T2|) . the points r„ should be all different modulo the symmetry. It suffices 
to choose them all on one sub-period of the base curve. 

Numerical exemple with the scalar field. The simplest base curve is 
a circle or radius a which is de-centered by ea with respect to the 2;-axis. 
We took e = 0.4, q = I and imposed the phase velocity Wph = 0.99999. The 
lowest modes are at ua = 5.0 and 7.0. They are stable at the 0.1 precision 
when C contains at least the set [—1,-1-5]. Fig.2 displays the cross section 
\I'(X, y) of the field ai z = t = 0, for ua = 5.0 mode. Note the symmetry 
f lT3|) . The maximum of |Re\E'| is not at the centre of the base, but farther 
from the helix axis. This is a manifestation of the centrifugal effect. 
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Figure 2: 3- dimensional representation of the scalar field \l'(x,?/) = 
$(0, x, y, 0) of the ua = 5.0 mode, for a circular base of radius a, off-centred 
by e = 0.4 from the axis, and q = 1. Left: real part; right: imaginary part. 

3 Propagating modes of the Maxwell field 

We gather the fields E and B in a 6-component vector F = {E,B}. In 
analogy with we look for solutions of the form 

F{t, X, y, z) = e^^^-^-* n[qz] T{X, Y) (14) 

with J'(X,y) = {£{X,Y),B{X,Y)}. n[qz] is the rotation matrix of ^ 
applied to the vectors S and B. The physical field is the real part of F. 



Boundary conditions. We introduce the vectors 

a(s) = r;,(s) = h/[(s) , f(s) = z + gz X rb(s) = I gxfe(s) I, (15) 





which are tangent to the wall at z = 0, and the normal vector u{s) = 
a{s) X r(s). The boundary conditions are 

{£r (Tr + S^a^ = 
qr£^ + £, = (16) 
cr^ + (qrB^ - B^)ar = 

The right set is written in cylindrical coordinates. These 3 conditions are 
not independent. Only two suffice. 
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Partial wave expansion. F can be expanded in free partial waves like $ 
in (IH]). However, since the photon has two polarisation states, there are two 
waves for each /, which now denotes the total (spin + orbital) angular momen- 
tum. We take the states of right and left chirality, or photon helicity x = =tl- 
The analogue of i^i,k{r) in is Jaj,p,«,x(^) = \^^uj,p,i,xi^),Buj,p,i,x{r)y For 
/ = 0, the cylindrical components are 

X^r) J 

(17) 

For / ^ 0, 




Suj,p,i,x = Si [ijjx-i - (p- x'(^fi^\+i] , ^u>,p,i,x = -'ix^uj,p,i,x- (18) 
V 2{p-x'uj)^x J 

Here si = sign(l), x' = siX: ^ = Kl and ip\ = J\{kr)/k^. These states have 
finite and linearly independent limits when — )■ cu^. The analogue of ([8]) is 

F(t,x,y,^) = e-^-* ^e^^^e^^W^ [ai :F^,pm,x=+sXr) + h :F^,pm,x=-sXr)\ -(19) 
I 

We made the choice to associate a/ and hi with positive and negative x' 
instead of x- Fo^' ^ = we set si = +1 in (fT^ . but we can alternatively take 
the (TM, TE) basis as in [7]. 



The point matching method in the Maxwell field case. Generalising 
the PWE-PM method of Section 2.2 we keep only values of / and write 
the first two conditions of f lT6|) for N points r„ of the base curve. It gives 2A'" 
linear equations linking N coefficients a/ and N coefficients bi. The 2A'" x 2A'" 
matrix of the system is given in [7] . The dispersion relation is the vanishing 
of its determinant. 



Numerical results for the Maxwell field. Taking the same de-centered 
circular base curve as in the scalar field case e = 0.4, but q = 1 and imposing 
Vph = 0.999, we obtain the lowest modes at ua = 2.28, 3.6, 4.8 and 5.0. They 
remain the same for various £ sets like [—6, +8], [—5, -|-9] and [—3, +11]. The 
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TM/TE character, measured by \Ez{0)/ Bz{0)\, are 0.3, 0.25 and 0.4. and 
5.4 respectively. The ua = 5.0 mode seems therefore more appropriate for 
electron acceleration. It has a figure of merit 

f.o.m. = |E^(0)|/[max(|E|) on the wall] ~ 1/8 (20) 

and a group velocity Vg = duj/dP = 0.68. Plots of the fields of this mode can 
be seen in [7j. Let us mention that \Br\ and look very similar to 
and \Sr\ respectively, which is reminiscent of a plane wave. 

Increasing the guide periods A = 27r/g while imposing fixed fph < 1 makes 
the u of each mode increase until the mode disappear at a; = oo. For q = 0.7, 
the two lowest modes are at ua = 2.7 and 4.2 with \Ez{0)/ Bz{0)\ = 0.25 and 
0.26 respectively. The lowest mode has f.o.m ~ 1/9. 

With a stronger helical bending of the waveguide (larger e and q), the 
mode frequencies decrease, but the convergence of the PWE-PM method 
becomes questionable. For e = 0.6, q = 1.5 and Vph = 0.999, three first 
modes are seen at coa = 2.2, 3.3 and 4.0 with \Ez{0)/ Bz{0)\ = 0.35, 0.3 
and 0.8 respectively. About twenty matching points were sufficient to obtain 
stable roots. However the field amplitudes were not stable. 

We also considered a non-circular base curve of the Pascal's limagon type 
(figure 1), parametrized by 

= a (1 — e cos s) , ipb = s — e' sins , s G [0, 27r] . (21) 

The inward bump is expected to give a stronger -^^(O). Indeed when a large 
part of the base curve is close to the origin, the boundary condition inhibits 
Ez{0). This is avoided with the bump. For e = e' = 0.3, q = 0.7 and 
Vph = 0.999, two modes are clearly seen at ua = 2.4 and 6.7. A state at 
coa = 3.5 was also seen for odd Igup, but it appeared only as a smooth dip of 
In |det{C„;}| for even kup- 

Finally we considered a base curve with two inside bumps and a rotational 
symmetry of order 2, paramerized by 

r;, = a (1 — ecoss) , (fb = {s — e' sms)/2 , sG[0,47r]. (22) 

The modes are either even or odd under this symmetry. Odd modes are 
not interesting for particle acceleration since Ez{0) = 0. Even modes have 
the property Ex{0) = Ey{0) = 0, which guaranties that electrons on the z- 
axis are not defiected and emit no synchrotron radiation. For e = e' = 0.1, 
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q = 0.7 and fph = 0.999, even modes are seen at ooa = 4.0 and 6.3, with the 
different / truncations [/inf,4up] = [~4,8], [—8,8] and [—10,10]. With more 
pronounced bumps {e = e' = 0.2 or 0.3), no stable root was found. 

4 Limits of the PWE-PM method 

We have encountered difficulties when the deformation of the wave guide 
relative to a straight cylindrical tube becomes too strong. This may be 
attributed not to the point matching method but to a lack of convergence of 
the partial wave expansion at r larger than a critical radius r^. 

Let us indeed consider the expansion ([HD, in the scalar field case, and 
suppose that it converges at some point of polar coordinates (ro,v5o)- Then 
it converges at all points (r, <f) with r < tq. It implies that the domain of 
convergence is bounded by a cylinder. 

Demonstration : The convergence of ([8]) at (ro, v^o) implies that \ai ipixi){''^o)\ 
is bounded by some number M. Using ([7]), we have 

\ai^pi,k(i)ir)\ < M 

for any r. At large we have ni \l\q — SiP with Si = sign(/) and the 
asymptotic form of the modified Bessel functions gives I\i\{Kir) / I\i\{KirQ) ~ 
exp[(|/|g — SiP) (r — Tq)]. Then 

|a/^/,fc(0(r)| <Mexp[s;P(ro-r)]el'l''(^-'-°) at large / . (24) 

If r < ro, the expansion ([H]) is bound in modulus by a convergent geometric 
series in |/|, therefore converges at I — )■ ±oo. 

If the base curve has a singularity, e.g., an anglar point of polar coordinate 
(ri,93i), the field is non-analytic at this point, therefore Vc < ri. On the 
other hand, for a strong enough deformation relative to the straight cylinder, 
part of the wave guide may be outside the cylinder of convergence, even 
if the base curve is analytic. A similar situation occurs when one uses the 
Rayleigh expansion for the wave diffracted by a grating : this expansion does 
not converge down to the bottom of the grooves if these ones are too deep 



(23) 
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5 Conclusion 



We have seen that an increase in the strength of the hehcal bending (increase 
of e and q in the case of a de-centered circular basis) leads to a decrease of 
the dimensionless frequencies ua of the low modes at imposed phase velocity 
^^ph < c (a characterizes the transverse size of the guide). One can also 
expect that it leads to a larger accelerating field -^2(0), at imposed power of 
the wave or imposed maximum field on the wall. We have indeed, in order 
of magnitude 

EJEt-^Pt/Pz, pT^l/a, Pz^i^ EJ Et l/{uja) . (25) 

Inward bumps of the base curve also favor large -^^(O). On the other hand, 
the PWE-PM method is limited to weak helical bending and sufficiently 
"round" base curves. Here we propose two numerical methods to find the 
propagating modes for larger bending and base curves with more pronounced 
structures. 



The finite element method. This is the standard way to solve numeri- 
cally the Maxwell equations. In the rotating frame these equations take the 
2-dimensional form at fixed u and P : 

iujBx = +dY£z - iP^Y - q^x + qiXdy - Ydx)£Y , 

iujBy = —dxSz + iP^x — Q^Y — liXdy — Ydx)£x , 

iuBz = +dx£Y - dySx , 

-iPBz = dxBx + dYBY + q{Ydx-XdY)B,, (26) 

and analogous equations obtained with the duahty transformation B ^ S, 
£ -B. 



Method of virtual currents. Instead of expanding the field in cylindrical 
partial waves, one does it with eigenfunctions of Pheiic (equation [5])which are 
regular inside the waveguide but singular at some outside places. A candidate 
is the retarded or advanced field produced by the charge- current distribution 



3^{t, v,z)=\]t\= e^'^^-^"* <5(') (r - n[qz]d) qv^,, z A r , (27) 
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with fph = uj/ P. TZ[qz] is the rotation matrix of ([T]). The x and y components 
of any vector are gathered in a bi-vector. The point {x,y,z) = (7l[qz]d,z) 
draws an hehx of path A = 27r/g intersecting the z = plane at point 
(d,0). p{t,x,y,z) is a periodic distribution of charges running along this 
helix with the longitudinal velocity Vz = fph and j is the associated current. 

is periodic in time with frequency u and is an eigenfunction of PhcUc 
with eigenvalue P, therefore the retarded or advanced field has the same 
property. Let us denote, for instance, by F*^*^^ = {E^'^^B^'^^} the retarded 
field. Different d's give different independent solutions. To take into account 
the fact that there are two polarisations, we also consider the dual field 
F^'^) = {B*^'^^ — E'^'^)}, which would be produced by a distribution of magnetic 
charge and current. In the z = plane these fields have a singularity in 
l/|r — d| at r = d ; they are Green's functions of fj26|) . Choosing singular 
points {di, d2, ■ ■ • d^, } external to the base curve, we may approximate the 
field of a mode by 



The coefficients and should be choosen to satisfy the boundary condi- 
tions on matching points ri, ...rjv of the base curve. 

Near a sharp inward bump of the base curve the field is strong and, 
by analytic continuation, one probably meets a singularity not far behind 
the curve, analogous to a point-like electrical image in electrostatics. This 
singularity may be simulated by placing one or a few d„'s on it or near to 
it. It amounts to replace the continuous surface charge and currents on the 
bump by charges and currents on virtual helical wires. 

Whichever method will be used, it will be usefull to compare it with the 
PWE-MP method in the case of weak helical deformation, as a check. 
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